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A search for the pair-production of heavy leptons (A°,L^) predicted by the type-III seesaw 
theory formulated to explain the origin of small neutrino masses is presented. The decay channels 
A° — >■ [I = e, [i, t) and — >■ W^u {v = I'lj., Vt) are considered. The analysis is performed 

using the hnal state that contains two leptons (electrons or muons), two jets from a hadronically 
decaying W boson and large missing transverse momentum. The data used in the measurement 
correspond to an integrated luminosity of 20.3 fb“^ of pp collisions at ^ = 8 TeV collected by 
the ATLAS detector at the LHC. No evidence of heavy lepton pair-production is observed. Heavy 
leptons with masses below 325-540 GeV are excluded at the 95% confidence level, depending on the 
theoretical scenario considered. 

PACS numbers: 14.60.Hi,13.35.Hb 


INTRODUCTION 

Experiments show that neutrinos have much smaller 
masses than charged leptons (see Ref. [T], and refer¬ 
ences therein). While in the Standard Model (SM) the 
charged fermions acquire masses by coupling to the Higgs 
{H) boson, the neutrinos may become massive via new 
physics beyond the SM, e.g. via the introduction of Ma- 
jorana mass terms [5]. These masses could be small due 
to the seesaw mechanism [am], which relies on new mas¬ 
sive states that couple to a charged lepton and the Higgs 
field. Among different models for the seesaw mechanism, 
the type-III model H ig introduces at least two extra 
triplets of fermionic fields with zero hypercharge in the 
adjoint representation of SU(2 )l that generate neutrino 
masses and couple to gauge bosons. This model predicts 
new charged and neutral heavy leptons that could be pro¬ 
duced in proton-proton collisions at the Large Hadron 
Collider (LHC). 

A search by the CMS experiment |6] excluded the type- 
III seesaw heavy leptons with masses in the range of 100- 
210 GeV, depending on theoretical assumptions. A re¬ 
cent search by ATLAS |7] also sets complementary limits 
on heavy leptons using the three-lepton final state. Sim¬ 
ilar searches have also been done by L3 experiment |8] 
ruling out charged heavy leptons with masses below 100 
GeV. 

In this paper, a search for heavy leptons predicted by 
the type-III seesaw mechanism is presented. The search 
explores the mass region above 100 GeV. A minimal 
type-III seesaw model [S] is used to optimize the anal¬ 
ysis strategy and interpret the search results. The model 
introduces a triplet with one neutral and two oppositely- 
charged leptons denoted by and L^, respectively. 
The heavy leptons decay into a SM lepton and a W, Z or 
Higgs boson. The heavy leptons are assumed to be degen¬ 
erate in mass. This assumption does not affect the result 
because in the case of a small mass splitting due to radia¬ 
tive corrections, the decays within the heavy leptons are 


highly suppressed nuj. The dominant production mech¬ 
anism for type-III seesaw heavy leptons in pp collisions 
is pair-production through the weak coupling to the W 
boson propagator: pp W* —?> N^L^, and the largest 
branching fraction is the one with two W bosons in the 
final state: —>■ {£ = e,/x,r) and —>■ 

(v = The production cross-section does not 

depend on the mixing angles between the SM leptons and 
the new heavy lepton states Va,{a = e, r), which en¬ 
ter only in the expressions for the L and N decay widths. 
The fraction of L and N decays to lepton flavor a is pro¬ 
portional to ba = IHcp/dUeP + + |HrP). The lim¬ 

its obtained may be interpreted in terms of a range of 
mixing angles and Yukawa couplings |5], allowing tests 
of a range of models with different couplings to gauge 
bosons and cross-section predictions mm- In the type- 
III seesaw model considered here, a benchmark point is 
defined by setting 14 to zero, so that bg and are deter¬ 
mined only by the ratio VejVfj,, taken to be 0.87 based on 
the separately-allowed maximum values of 14 and I 4 in 
Refs. [ISHIl]. This choice results in values of b^ = 0.53, 
bfj, = 0.43 and br = 0. The search is performed for the 
process pp —?> N'^L^ —where one W bo¬ 
son decays leptonically and the other W boson decays 
hadronically, resulting in a lepton pair in the final state 
with either the same charge (same-sign, SS) or with the 
opposite charge (opposite-sign, OS). 


DATA SAMPLE AND MONTE CARLO 
SIMULATION 

The analysis uses data from ^/s = 8 TeV pp colli¬ 
sions at the LHC that were recorded by the ATLAS 
detector using single-electron and single-muon triggers. 
A detailed description of the ATLAS detector can be 
found elsewhere [ig. The data sample corresponds to 
20.3±0.6fb“^ [T7] of integrated luminosity. Data quality 
criteria are applied to ensure that events were recorded 
with stable beam conditions and with all relevant sub- 
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detector systems operational. The triggers [l8] are fully 
efficient for leptons with px > 25 GeV, where transverse 
momentum px is defined as the magnitude of the momen¬ 
tum component orthogonal to the beam axis. Events are 
required to have a reconstructed collision vertex with at 
least three associated tracks, each with px > 400 MeV. 
In events with multiple vertices, the vertex with the 
largest ^ Px of associated tracks is taken as the primary 
event vertex. 

Monte Carlo (MC) samples are used to optimize the 
event selection and to model the kinematics and nor¬ 
malization of most background processes. Signal sam¬ 
ples are generated for heavy lepton masses in the range 
100-600 GeV. MadGraph 5 [12] is used to calcu¬ 
late the matrix elements for each process, while MadE- 
VENT [3D] with the MSTW2008 parton distribution func¬ 
tions (PDF) set |2T] simulates the initial hard scatter¬ 
ing and the N and L decay. Pythia 8.153 [22] is used 
to simulate the decays of W bosons and the underlying 
physics by providing parton showers and hadronization, 
as well as adding initial- and final-state radiation (ISR 
and FSR) to the events simulated in MadEvent. The 
main background sources arise from the production of 
a Z boson in association with jets (Z-|-jets), single and 
pair-production of top quarks, and diboson production 
(WWjWZ, ZZ). The Z-|-jets and diboson processes are 
simulated with Sherpa 1.4.1 [3Dj, a generator based on a 
multileg matrix element calculation matched to the par- 
ton shower using the CKKW prescription |33| , and using 
the CTIO [3^ PDF set. For diboson production, both 
the electroweak and strong production processes are sim¬ 
ulated [3^. Top-quark pair events and single-top-quark 
events in the VFt-channel and s-channel are simulated 
using MC@NLO 4.03 [27], which is interfaced to Her- 
wiG 6.520 [35] and Jimmy 4.31 [3D] with the CTIO PDF 
set. Top pair production in association with a vector bo¬ 
son, tt + WjZ, is simulated using MadGraph with the 
CTEQ6L1 PDF set [30], interfaced to Pythia 8.153 for 
parton showering and hadronisation. Single top quark 
production in the t-channel is simulated using AcerMC 
v3.8 [S] with Pythia 6.426 [33] and the CTEQ6L1 PDF 
set. 

All samples of simulated events include the effect 
of multiple pp interactions in the same and neighbor¬ 
ing bunch crossings (pileup) by overlaying simulated 
minimum-bias events on each generated signal and back¬ 
ground event. The number of overlaid events is chosen to 
match the average number of interactions per pp bunch 
crossing observed in the data as it evolved throughout the 
data-taking period (giving an average of 21 interactions 
per crossing for the whole data-taking period). The gen¬ 
erated samples are processed through the GEANT4-based 
detector simulation [33] [53] or a fast simulation using a 
parameterization of the performance of the calorimetry 
and Geant4 for the other parts of the detector [35] . The 
standard ATLAS reconstruction software is used for both 


simulated and collision data. 


OBJECT DEFINITIONS 

The reconstructed objects used in this analysis are 
electrons, muons, jets, and missing transverse momen¬ 
tum. Electrons are reconstructed from clusters of energy 
depositions in the calorimeter that match a track recon¬ 
structed in the inner detector (ID) and satisfy the “tight” 
criteria defined in Ref. [36] . The electrons are required to 
havepx > 25 GeV and pseudorapidity |?7| < 2.47 [37], ex¬ 
cluding the transition region between the barrel and end- 
caps in the liquid argon calorimeter (1.37 < \rj\ < 1.52). 
Muons are reconstructed by combining ID and and muon 
spectrometer tracks that are spatially matched and have 
consistent curvatures. The muon tracks are required to 
havepx > 25 GeV and \r]\ < 2.5. In addition, leptons are 
required to be isolated from other tracks and calorimet¬ 
ric activity [35]. To ensure that leptons originate from 
the interaction point, requirements of \dQ\/ado < 3 and 
l^osin^l < 0.5 mm are imposed on the electrons and 
muons, where do (-^o) is the transverse (longitudinal) im¬ 
pact parameter of the lepton, and cr^o is the uncertainty 
on the measured do- The lepton impact parameters are 
measured with respect to the event primary vertex. 

Jets are reconstructed from three-dimensional topolog¬ 
ical clusters of energy depositions in the calorimeter us¬ 
ing the anti-fct algorithm [39] with a radius parameter 
of i? = 0.4. The energies of jets are calibrated to the 
hadronic energy scale by correcting for energy losses in 
passive material, the non-compensating response of the 
calorimeter, and extra energy due to multiple pp inter¬ 
actions [40]. The jets are required to have px > 30 GeV 
and I77I < 2.8. For jets with px < 50 GeV and \r]\ < 2.4, 
the summed scalar px of associated tracks from the re¬ 
constructed primary vertex is required to be at least 25% 
of the summed scalar px of all tracks associated with the 
jet. In the pseudorapidity range \r]\ < 2.5, jets containing 
&-hadrons are identified using a 6-tagging algorithm [41] 
with an efficiency of 70% and with a misidentification 
rate for selecting light-quark or gluon jets of less than 
1%. The identification efficiency of the algorithm for 
jets containing c-hadrons is 20%. The efficiencies and 
misidentification rates are determined from ti MG events. 

The missing transverse momentum vector (with its 
magnitude Ax”***) is derived using the calorimeter cell 
energies within \r]\ < 4.9 and corrected on the basis of 
dedicated calibrations of the associated physics objects 
including muons [42] . Calorimeter cells containing energy 
depositions above noise and not associated with high-px 
physics objects are also included. 
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EVENT SELECTION 

Events that contain exactly two reconstructed leptons 
(electrons or muons), at least two jets, and no 6-tagged 
jets are selected. One of these leptons is required to 
match the object upon which the event was triggered. 
Different sets of optimized selection criteria are used for 
the events in the OS and SS final states. The optimiza¬ 
tion is done using simulated heavy lepton pair-production 
events at a benchmark mass of 300 GeV. For the OS (SS) 
final state, the leading and next-to-leading lepton candi¬ 
dates are required to have px greater than 100 (70) GeV 
and 25 (40) GeV, respectively. The invariant mass of the 
two lepton candidates is required to be larger than 130 
(90) GeV in order to suppress background from the pro¬ 
duction of Z-l-jets. The hadronically decaying W can¬ 
didate is formed by combining the two jets with high¬ 
est pt, and the pt of the first and second leading jets 
are respectively required to be larger than 60 (40) and 
30 (25) GeV, for the OS (SS) final state. The invari¬ 
ant mass of the W candidate, rnjj, is required to be 
between 60 and 100 GeV. Events selected in the OS 
(SS) final state are required to have a of at least 

110 (100) GeV and, for OS events, an angular separation 
A-Rjj = \l (A? 7)2 -t (A()))2 < 2 between the two jets with 
highest pt: where Ap and At)) are defined as the differ¬ 
ences in pseudorapidity and azimuthal angle between the 
jets. 

BACKGROUND ESTIMATE 

The background in this search can be classified into 
two categories based on the origin of the charged lepton 
candidates. The first category of backgrounds consists 
of events in which two leptons are produced via the de¬ 
cays of W or Z bosons and are correctly reconstructed. 
This category of backgrounds includes the production 
of Z-l-jets, ti, Wt single-top-quark and diboson events. 
Smaller contributions originate from ti + WjZ events. 
Gontributions from triboson events, such as WWW^ and 
events containing a Higgs boson, are negligible. The 
number of events from this background category is esti¬ 
mated using the simulated samples described previously. 

The second category corresponds to all other sources, 
such as events containing at least one particle that is 
incorrectly identified as a lepton, or a lepton which origi¬ 
nates from secondary interactions and decays, which are 
together denoted as fakes, and events with a lepton whose 
charge is incorrectly determined. Fake electrons originate 
primarily from jets that have a large electromagnetic en¬ 
ergy fraction passing the electron selection requirements, 
photon conversions, and electrons from semileptonic de¬ 
cays of charm or bottom hadrons. Fake muons include 
muons arising from semileptonic decays of charm or bot¬ 
tom hadrons, in-flight decays of pions or kaons, or ener¬ 


getic particles that reach the muon spectrometer. 

For the OS final state, the background contribution is 
dominated by events from the first category. For the SS 
final state, the expected background is very small and 
is also primarily from the first category. Events from 
Z-|-jets production can mimic the SS signal if the charge 
of one of the electrons is incorrectly determined. The 
background events from this contribution are modeled 
by simulation, with correction factors derived using Z —>■ 
ee events in data. The probability of misidentifying the 
muon charge is negligible. 

The background contribution in the second category 
(i.e. fake leptons) is estimated using an in situ tech¬ 
nique [38] that relies minimally on simulation. This is 
done by reweighting a complementary set of events, se¬ 
lected by changing the electron identification criterion 
from “tight” to “loose” |13] and by loosening the muon 
Mo|/<7do and the electron and muon isolation require¬ 
ments, while keeping the event selection otherwise iden¬ 
tical. The reweighting factors are defined as the ratio of 
the number of events containing a lepton that satisfies 
the nominal criteria to the number of events containing 
a lepton that only fulfill the relaxed criteria. These fac¬ 
tors are measured as a function of the candidate pt and 
r] in data samples that are enriched in fake leptons [38] . 
Gorrections to the factors due to true leptons from vec¬ 
tor boson decay in the background-enriched samples are 
taken from MG simulation. 

Figure [2 shows a comparison of the missing transverse 
momentum distribution of data, expected backgrounds 
and signal predictions when all the selection require¬ 
ments, except for the missing transverse momentum re¬ 
quirement, are applied. The shape and the rate of the 
background estimate is in good agreement with the data. 

The background estimates are validated by comparing 
the predicted numbers of events in simulation to those 
observed in data in several control regions that have event 
selection criteria similar to those for the signal region. 

The control region for top quark pairs is defined by 


TABLE I. Event yields for opposite-sign (OS) and same-sign 
(SS) selection for predicted backgrounds, data, and type-III 
seesaw lepton pair-production with masses of 150 and 300 
GeV. The reported errors include both the statistical and sys¬ 
tematic uncertainties. 




OS 



SS 


Fake Leptons 


1.4 ± 

0.9 

0.67 

± 

0.42 

Z-|-jets 


2.4 ± 

1.2 

0.06 

± 

0.23 

ww/wz/zz 


9.2 ± 

2.9 

1.95 

± 

0.58 

tt i+W/Z) and 

single top 

1 17.9 ± 

6.9 

0.47 

± 

0.25 

Total 


31.0 ± 

7.7 

3.15 

± 

0.80 

Data 


25 
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Signal m^/jv = 

150 GeV 

9.5 ± 

1.6 

20.3 

I i 

2.3 

Signal rtiL/N = 

300 GeV 

12.2 ± 

0.6 

5.7 

± 

0.5 
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FIG. 1. The missing transverse momentum distribution of opposite-sign (left) and same-sign (right) events for data and 
predictions. The dotted and dashed lines show the expected distribution of type-III seesaw lepton pair-production with masses 
of 150 and 300 GeV, respectively. The events shown here are required to pass all selection requirements except that on the 
missing transverse momentum. The uncertainties shown include both statistical and systematic uncertainties. 


selecting events with two 6-tagged jets. In this region, 



OS ee SS ee OS pp SS pp OS ep SS ep 


FIG. 2. Event yields for opposite-sign (OS) and same-sign 
(SS) selection for ee, fifj,, and e/r predicted backgrounds, 
data, and signal events featuring type-III seesaw lepton pair- 
production with masses 150 and 300 GeV. The reported un¬ 
certainties include both the statistical and systematic uncer¬ 
tainties. 


according to MC simulation, all the events are from top 
quark pair and single top quark production with a neg¬ 
ligible contribution from other sources. MC simulation 
predicts 26 ± 3 (stat.) events and in data 32 events 
are observed. The scale factor, the ratio of the observed 
and predicted event yields, is found to be consistent with 
unity. 

The diboson control region is obtained using a WZ- 
enriched sample of events containing three leptons with¬ 
out a requirement on the missing transverse momentum. 
In this control region, according to MC simulation all 
the events are from diboson production with a negligible 
contribution from other sources. MC simulation predicts 
II ± I events and in data 9 events are observed. The 
scale factor is found to be consistent with unity. 


SYSTEMATIC UNCERTAINTIES 

The uncertainties on the rate of top quark and diboson 
backgrounds due to potential differences between data 
and MC simulation are evaluated using the statistical 
uncertainties of the measured scale factors in the control 
samples. They are the dominant systematic uncertain¬ 
ties (^ 35%) on the background estimates. For top quark 
production, an additional systematic uncertainty is con¬ 
sidered to account for a potential difference between scale 
factors in the control region and the signal region. This 
is done by comparing the nominal tt sample to alterna- 
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FIG. 3. Left: Observed (solid line) and expected (dashed line) 95% C.L. upper limits on the type-III seesaw heavy lepton cross- 
section as a function of the heavy lepton L or N mass assuming 6e = 0.43, 6^ = 0.57, and br = 0. The bands surrounding the 
expected limit correspond to one and two standard deviations on the expected limit. The large-dashed (dot-dashed) line shows 
the theoretical prediction for m/,/jv-dependent (maximal) branching fraction for decays to a IT boson. Right: The expected 
(large-dashed line) and observed (shaded region) upper limits at the 95% C.L. on the 577(1,^ —>■ Wu) x BR{N‘^ —>■ VF^/^) 
versus The dashed line corresponds to one standard deviation around the expected limit. The dotted line shows the 

nominal mass-dependent branching ratio of the type-III seesaw model. 




FIG. 4. Observed (solid line) and expected (dashed line) 95% C.L. upper limits on the type-III seesaw heavy lepton cross-section 
as a function of the heavy lepton L or N mass for exclusive coupling to electrons (left) and exclusive coupling to muons (right). 
In the left plot, the dashed line is covered by the solid line. The large-dashed (dot-dashed) line shows the theoretical prediction 
for dependent (maximal) branching fraction for decays to a IF boson. 


five ti MC samples. These samples include events that 
are generated using Powheg-Box 1.0 (patch 4) |44Il46] 
and the leading-order, multileg generator Alpgen v2.13 
m- The Powheg-Box generator is interfaced to the 
Pythia 6.426 showering routines with either CTIO or 
HeraPDF [48] PDF sets and with the Powheg /idamp 
parameter set to either the mass of the top quark or in¬ 
finity [35]. Alpgen is interfaced to Herwig 6.520, and 
used to simulate top pair events with up to four addi¬ 
tional partons in the matrix element. The uncertain¬ 
ties due to QCD ISR and FSR modeling are estimated 
with samples generated with AcerMC v3.8 interfaced 


to Pythia 6.426 in which the parton shower parameters 
are varied in a range consistent with a measurement of 
additional hadronic activity in ti events jSnj. The differ¬ 
ences observed in the signal region by using different MC 
simulations are about 35%. 

For the Z+n-jets (n > 2) background estimate, the 
dominant systematic uncertainty in the OS final state 
is from the uncertainty on its production cross-section 
(^ 50%) |5T]. For the SS final state, the systematic un¬ 
certainty is dominated by the statistical uncertainty on 
the measured electron charge misidentification rate. 

Uncertainties on the background estimate due to fake 
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leptons are determined in dedicated studies using a com¬ 
bination of simulation and data. They account for poten¬ 
tial biases in the method used to extract the reweighting 
factors, and for the dependency of the reweighting factors 
on the event topology. 

For both the predicted signal and background event 
yields, uncertainties resulting from detector effects from 
jet energy scale and resolution [40) . lepton reconstruction 
and identification efficiencies [361152] , lepton momentum 
scales and resolutions |5ll [S3], and missing transverse 
momentum [42] are considered. They are typically small 
(1-5%). The theoretical uncertainties on the signal pro¬ 
duction cross-section and acceptance, such as PDF choice 
and ISR and FSR modeling, are found to be negligible. 

The background estimates and their uncertainties are 
tested in two other regions: a Z-|-jets control region and 
the hadronic W sidebands. The .Z-l-jets control region is 
selected by requiring the invariant mass of the oppositely- 
charged lepton pair to be consistent with the Z boson 
mass. In this region, approximately 70% of the events 
are from a Z boson produced in association with jets. 
The predicted number of events is 34 ± 4, where the 
error includes both the statistical and systematic uncer¬ 
tainties, and 32 events are observed in data. For the 
hadronic W sidebands, the analysis was repeated using 
the same event selection, but requiring the invariant mass 
of the hadronically decaying W candidate to be 35 < 
rrijj < 60 GeV or 100 < mjj < 125 GeV. This selection 
provides samples dominated by background events with 
kinematic properties similar to those of the signal candi¬ 
dates. In this region the predicted number of events is 34 
± 7 events, where the error includes both the statistical 
and systematic uncertainties, and in data 18 events are 
observed. The data are in agreement with the predictions 
within 1.9 standard deviations. 


RESULTS AND INTERPRETATION 

Table |l| shows the predicted numbers of signal and 
background events and the observed data events in the 
signal region. The data agree with the background-only 
hypothesis. Figure shows the flavor composition of 
the simulated signal and background events, and of the 
observed events in data. In the absence of any signif¬ 
icant data excess, upper limits on the production rate 
of pp —>■ —)■ at the 95% confidence 

level (G.L.) are derived as a function of the heavy lepton 
mass using the GLs method [S3]. The results of limit 
calculations combining the observations in the OS and 
SS final states are shown in Figs. and [^ With the de¬ 
fault mixing angles considered here, heavy leptons with 
masses less than 335 GeV are excluded by the analysis, 
while masses less than 475 GeV are excluded in the sce¬ 
nario in which heavy leptons can only decay to the W(. 
ov Wy final states. For comparison, upper limits are also 


calculated for different theoretical assumptions, such as 
exclusive coupling between the heavy leptons and muons 
{he = 0,5^ = 1) or electrons {be = 1,6^ = 0 ). In the 
limit calculations for exclusive couplings, events in the 
two-muon (two-electron) final state are excluded for ex¬ 
clusive electron (muon) coupling. Heavy leptons with a 
mass below 400 (325) GeV can be ruled out by the data 
in the case of exclusive coupling to muons (electrons). 
Masses less than 540 (470) GeV are excluded in the sce¬ 
nario in which heavy leptons can only decay to the W£ 
OT Wy final state. 


CONCLUSIONS 

A search for the pair-production of heavy leptons pre¬ 
dicted by the type-III seesaw model is presented. The 
analysis is performed using a final state that contains two 
leptons, two jets from a hadronically decaying W boson, 
and large missing transverse momentum. The data used 
in the search correspond to an integrated luminosity of 
20.3 fb“^ of pp collisions at ,/s = 8 TeV collected by the 
ATLAS detector at the LHC. No evidence of heavy lep¬ 
ton production is observed. Heavy leptons with masses 
below 325-540 GeV are excluded at the 95% confidence 
level, depending on the considered theoretical scenarios. 
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